Scaling properties of g-breathers in nonlinear acoustic lattices 
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Recently g-breathers - time-periodic solutions which localize in the space of normal modes and 
maximize the energy density for some mode vector go - were obtained for finite nonlinear lattices. We 
scale these solutions together with the size of the system to arbitrarily large lattices. We generalize 
previously obtained analytical estimates of the localization length of g-breathers. The first finding is 
that the degree of localization depends only on intensive quantities and is size independent. Secondly 
a critical wave vector k m is identified, which depends on one effective nonlinearity parameter, q- 
breathers minimize the localization length at fco = k m and completely delocalize in the limit ko — > 0. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 
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INTRODUCTION 

Spatially extended nonlinear Hamiltonian systems 
serve as starting models for the study of excitations in 
many branches in physics, e.g. anharmonic vibrations 
of crystal lattices, mesoscopic and nanoscopic systems, 
molecules, but also of electromagnetic, acoustic and other 
waves in nonlinear media, to name a few. They have 
been studied over many decades in order to understand 
such intriguing material properties as heat conductivity, 
thermal expansion, turbulence, confinement of light, but 
also general mathematical aspects such as thermaliza- 
tion, mode-mode interactions, etc. While in any realis- 
tic situation damping and energy input have to be con- 
sidered as well, these dissipative effects are often weak 
enough to allow the observation of the underlying Hamil- 
tonian excitations. 

Recently it was shown pj, that a one-dimensional an- 
harmonic atomic chain allows for exact time-periodic so- 
lutions which localize exponentially in the space of nor- 
mal modes and have their maximum energy on a mode 
with mode number go . A g-breather, being periodic in 
time, can be viewed as one normal mode which is dressed 
by several other normal modes in a neighbourhood of go 
and has an infinite lifetime. The existence and properties 
of these g-breathers allowed to explain all major ingre- 
dients of the famous Fermi-Pasta-Ulam problem (FPU) 
0, 01, 0| in a clear and constructive way. The FPU prob- 
lem concerns the nonequipartition of normal mode ener- 
gies on time scales which can be many orders of mag- 
nitude larger than the characteristic vibration periods. 
The key ingredient for the construction of g-breathers 
is a finite nonlinear system Q]. This is a very general 
condition and may apply to many other systems as well. 
It has been recently successfully tested by considering 
FPU models with lattice dimension d — 2,3 0] and also 
discrete nonlinear Schrodinger models (DNLS) in vari- 
ous lattice dimensions Q. Further studies also revealed 
the persistence of g-breathers in thermal equilibrium pj, 



which shows their relevance for statistical properties of 
extended systems. Previous studies of the FPU problem 
suggested that the effect of nonequipartition will disap- 
pear for large system sizes 0, Q . That seems to imply 
a disappearance of g-breathers in that limit. Here we 
show that g-breathers persist and have invariant proper- 
ties for large system sizes. These results are particularly 
important because they apply to macroscopic systems. 

Consider a generic model of a D-dimensional nonlinear 
lattice of size N\ X • • • X Njj , defined by a Hamiltonian 
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where x n , p n are canonical variables. U(x) and V(x) 
are anharmonic on-site and interaction potentials, respec- 
tively. Their Taylor expansion around x = starts with 
quadratic terms, n = {m, . . . , Tin} is a £>-dimensional 
lattice vector with n; = 1, JVj. e; denotes a unitary lat- 
tice vector along the dimension I. Note that x n = is an 
equilibrium state of the system. We will consider the case 
of fixed boundary conditions (BC). We have also studied 
free and periodic BC with similar results. 

This Hamiltonian can be expressed in terms of nor- 
mal modes P q , Q q of the linearized problem, obtained 
by skipping all anharmonic terms in the potentials: 



- (P 2 



(2) 



where E q is the energy of a given normal mode and H mt 
is the mode interaction part of the Hamiltonian which 
appears due to anharmonicity. The integer components 
of the mode vector q enumerate the modes. That is a 
class of models for which exact g-breather solutions may 
exist Q, 0|. Such solutions are time-periodic and the 
normal mode energies are exponentially localized around 
a mode vector q . 

We search for a way of scaling a given solution Q q (t) 
of a system Q of size N to a solution Q q (t) of a system 
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FIG. 1: Constructing a solution in a chain of double size 
in real space. Boundaries and the additional lattice site in 
the center are marked with crosses. Initial size is N=32. The 
momentary displacements versus lattice site are shown, (color 
online) 

of scaled size. It consists of scaling the values of mode 
variables ^/r times and scaling their indices r times, filling 
the gaps with zeros: 

where we omitted all mode indices except the Z-th com- 
ponent. For fixed BC qi — l,Ni, and the scaled system 
size is assumed Ni + 1 = r(Ni + 1). The phase space of 
the scaled system then possesses an invariant subspace 
(also coined a bush of modes 0|). This procedure can be 
applied to construct a solution to a system of arbitrarily 
large size, increasing r. In order to identify the applica- 
bility of (jHJ, we will use a renormalization in real space 
which is equivalent to ©. 

Consider D = 1 in Q . The equations of motion read 

Xn = f{x n ) + P(x n +1 ~ X n ) - lf(x n ~ X n -l) ■ (4) 

Here f(x) = -U'(x), tp(x) = V'(x) and n = TJST. For 
fixed BC xq = xjv+i = 0. We further assume an odd 
restoring force f{x) = —f{—x). Then <@J is invariant 
under the following combined parity and sign reversal 
symmetry: x n — > — xn+i-u- The transformation Q is 
given by an alternation of spatial blocks, obtained from 
the previous by parity and sign reversal transform. The 
blocks are separated by additional nonexcited lattice sites 
(see Fig.Q): 

X n (t) = {x\(t), . . . ,XN(t),0, -XN(t), . . . , ~Xl(t), 0, . . . }. 

It is straightforward to observe that if x n (t) is a solu- 
tion to the initial system, then x n (t) is a solution to the 
scaled-size system. The scaling rule <EJ is thus confirmed 
for ID chains Q. 

We generalize the above results to higher lattice dimen- 
sions. The transformation to mode variables is a super- 
position of ID transforms. Then, the transformation <EJ 
along a lattice direction I corresponds to the real-space 
transform already discussed for the D — 1 case. It yields 
a solution to a scaled-size system if both f(x) and (p{x) 
are odd functions. 

Given a g-breather solution for the original finite sys- 
tem, we can thus scale the g-breather solution to larger 



system sizes. Its total energy is scaled in all cases like 
E = rE, which is ensured by the block structure of 
the scaling and the local structure of the coupling in 
the Hamiltonian JQ. The time-dependent mode ener- 
gies E q are transformed as E q (t) = rE q (t) for q = rq, 
and E q — for other q. Introducing the energy density 
as e = E/(N + 1) and wave number k = nq/(N + 1), it 
is straightforward to observe that the scaling procedure 
leaves the energy density and the wave number (vector) of 
a g-breather invariant. Together with the rigorous proof 
of existence of g-breathers for finite systems [l( we arrive 
at a rigorous proof of existence of these excitations for in- 
finite system sizes, with proper scaling and under certain 
restrictions for the potential functions. Since the scaled 
solutions are embedded on mode bushes Q, the question 
arises whether q-breathers with other (or any) values of 
ko exist in large systems as well. The fact that the scaling 
preserves the localization properties of the scaled excita- 
tions suggests a positive answer. Below we will address 
this question. We also note that we needed certain sym- 
metry properties of the potentials for the scaling to work. 
It however does not imply that for cases with less sym- 
metries (/-breathers do not exist. 

The localization properties of q-breathers in the nor- 
mal mode space have been obtained analytically using 
asymptotic expansions for FPU-models in various lattice 
dimensions. For / = and ip = x + fix 3 , the result from 
1, 5| is expressed in total energies and mode numbers 
£(2n+i) 90 = A"S 90 with V\ = 3f3Eg (N + l)/(8n 2 q 2 ) for 
D = 1 and qo <^ N. We conclude that it must be pos- 
sible to substitute densities and wave numbers instead, 
and obtain an expression which is independent of the ac- 
tual system size. Indeed the outcome of our substitution 
is written in the following way: 

lne fc = ^-l)lnVA + ln £fc0 , ^A = ^^. (5) 

Note that these results hold actually for D — 1,2,3 jjj. 
We observe the same scaling properties of © as sug- 
gested by the real space renormalization procedure from 
above. Equation © implies a smooth dependence of A 
on Ek and ko. We thus test whether the scaling can 
be observed, when the size of the system takes differ- 
ent values, and the new wave number ko is chosen to be 
the nearest one to the original ko value. We compute 
g-breathers for that model with D = 1 for various sys- 
tem sizes starting with 7V = 15 and plot the numerically 
found A as a function of N in the left panel in FigEJ The 
mode number q is chosen to be the nearest integer to 
ko(N + l)/ir, ko = 7r/15. We obtain A from the ratio of 
the energy densities E5 qo / e^ qo . First we observe that A is 
independent on N when N + 1 = r(No + 1). Secondly we 
observe fluctuations of A around a mean value for other 
values of N due to the fact that for these system sizes 
the closest wave number to ko will nevertheless be slightly 
different. These deviations decrease with increase of the 
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FIG. 2: Dependence of the localization parameter A on the 
system size N of a q-breather (left panel) for the /3-FPU chain 
at constant energy density, e = 4- 10~ 4 and 2- 10~ 3 (for lower 
and upper curves, respectively) and /3 = 1. go is chosen to be 
the nearest integer to ko{N + 1)/ty, fco = n/15. Open circles 
mark system sizes N + 1 = r(No + 1), No = 15, for which the 
prediction is that A should be invariant on the integer r. The 
right panel shows the smooth dependence of the measured A 
values on the actual wave numbers fco used in the left panel, 
(color online) 



system size, and thus the fluctuation amplitude in FigO 
decreases as well. The piecewise smooth curves show 
that A depends smoothly on fco, since by construction we 
probe A with wavenumbers slightly varying around fco. 
That is nicely confirmed in the right panel in FigEl 

Let us analyze the fc-dependence of 101 for D = 1. For 
g-breathers at fixed average energy density e it follows 
within the approximation of exponential localization that 
£ k — (1 — ^) £ - Together with the definition of A in © 
we find, that the inverse localization length in fc-space is 
given by the absolute value of the slope S : 



3/3 



fco 2v £ /Kq 8 

(6) 

A depends both on fco and Ek . S vanishes for fco — > and 
has its largest absolute value max(\S\) « 0.7432/^ at 
kmin ~ 2.577ia For a fixed effective nonlinearity param- 
eter v the g-breather with fco = fc mm shows the strongest 
localization. Especially for fco — > the g-breather de- 
localizes completely. With increasing v, the localization 
length of the g-breather for fco = fc m m increases. For 
fco 3> v it follows S ~ 2/fco ln(z//fco) and for fco <C v we 
find S « -fc /(2;y 2 ). 

We plot in FigEl the dependence S(k ) for the /3-FPU 
chain at three different energy densities. If v or N are 
small enough, then the first nonzero fco value will appear 
for fco > fcmm. Increasing v or N we shift some allowed 
low lying fc values to the left of the minimum fco < fc m m • 
For very large systems (dense filling of the x-axis in FigEl 
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FIG. 3: The slope S as a function of k for the /3-FPU 
chain and three different energy densities £ = 6.08- 10" 4 , 9.6- 
1(T 4 ,1.57 ■ 1(T 3 and f3 = 1 (dashed curves, from bottom 
to top). Symbols and eye-guiding solid lines: estimate of the 
slope from numerical computations of g-breathers for iV = 149 
and iV = 359. (color online) 
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FIG. 4: The master slope function S m (z) (dashed line). Dif- 
ferent symbols correspond to the numerically obtained slopes 
from FigEl which are scaled accordingly and fall on a single 
curve, (color online) 



with allowed fco values) we thus expect that among them 
there is an optimum wavelength which provides strongest 
localization. We test our prediction by computing the 
slope for various g-breathers of the /3-FPU chain. The re- 
sults are shown in FigEl (symbols) . We nicely observe an 
optimum value of fco for localization. Increasing the non- 
linearity parameter (by either increasing the nonlinearity 
strength or the energy density) the critical wavenumber 
is shifted further away from the edge of the spectrum and 
gets shallower, as predicted. Deviations from the theo- 
retical curves for small fco are due to strong nonlinear 
corrections to our estimates, while deviations at larger 
fco are finite size corrections to the analytical estimates. 
Note the smooth dependence of S on fco which does not 
depend on the system size, thus confirming the scaling 
results from above. We plot in FigQ] the single master 
slope function S m {z) = vS with z — fco/V and scale all 
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FIG. 5: The profiles of g-breather solutions for iV = 359, 
/3 = 1, e = 9.6 ■ 10 -4 and three different values of go which 
yield ko/ir = 0.011,0.022,0.067 and correspond to the most 
left point on the middle curve in Figl^J its minimum and a 
point to the right of it. Symbols denote the mode energies 
at the moment of vanishing coordinates xi = 0. The vertical 
bars denote the range of mode energy values taken during one 
period of the (/-breather, (color online) 

numerically obtained slopes as well. The numerical data 
all condense on one single curve at small fco . Thus higher 
order corrections to the decay law do not alter the scaling 
properties of g-breathers in the limit of small fco. 

The profiles of three g-breather solutions for s = 
9.6 • 10 -4 are shown in FigEJ The symbols correspond 
to normal mode energies at the moment when all coor- 
dinates vanish. Among them is the g-breather with the 
strongest localization. The precision of computation is 
10~ 8 (see P| for details). Since the normal mode energies 
are not conserved quantities, we show as well the fluctu- 
ation range for each of them. These fluctuations become 
stronger at particular fc- values and are possibly due to a 
closely nearby lying resonance, which nevertheless does 
not destroy the localization profile. 

Let us discuss the obtained results. The reason for the 
weaker localization of g-breathers when fc ^S> v is the in- 
creasing distance 3fco between modes excited in consecu- 
tive orders of perturbation theory. The delocalization for 
ko — > however is due to an approaching of resonances 
nu>k — > for some integer n Note that the same 
approaching of resonances holds at the upper frequency 
cutoff where the frequency detuning is quadratic in k and 
the relevant integer n = 1 . We computed the dependence 
of the slope S there and obtained a behaviour similar to 
the one in FigEI We expect the above results to qualita- 
tively hold independent of the dimension D. It remains a 
challenging task to perform computations for e.g. D = 2, 
since one needs about 10 5 lattice sites, which is presently 



not reachable with our numerical tools. 

We fixed the average energy density in order to ensure 
finite temperatures. If the energy density Ek Q is fixed, 
then g-breathers will delocalize for some nonzero value 
of ko. For a finite lattice e f» Nek at that point. 

We derived similar results for the ev-FPU model with 
(p = x + ax 2 using the estimates from pj. We obtain 
strongest localization for k m w 2.39£ and max(\S\) w 
1.5/| with the effective nonlinearity parameter £ = 

y/a/^e 1/4 . For small k the slope S (2fc /£) 1/3 . 

In all cases the localization becomes meaningless when 
| ^S' | 1 « 7r which is the size of the first Brilloin zone. 
That happens at fc ~ 3/3e/(47r) (/3-FPU) and at 
fc « a 2 e/(2n 5 ) (a-FPU)._Well defined and localized q- 
breathers exist for fco 3> fco- Strong resonances destroy 
them for fco <C fco and lead to an effective redistribution 
of mode energy in that hydrodynamic regime. The same 
reasoning defines a critical nonlinearity value, for which 
k m reaches the center of the band. It can be roughly esti- 
mated as v, £ ~ 1. For larger values of v, £ the system will 
enter a regime of strong nonlinearity, where g-breathers 
may become meaningless. 

We considered standing waves. We expect the results 
to be also of importance for travelling waves which are 
reflected at boundaries or inhomogeneities. 
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